Inverse modeling of lichen growth curves and implications for lichenometric dating     

《Quaternary Geochronology》2022

Lichenometric dating represents a quick and affordable surface exposure dating method that has been widely used to provide a minimum age constraint on tectonic and geomorphic landscape changes as well as buildings and anthropogenic landscape changes in various settings during the late Holocene. Despite its widespread usage, this method has several limitations. Major problems relate to the sampling of lichen population on any given rock surface and the modeling of growth curves. In order to overcome these issues, it has been suggested to subdivide the rock surface into some areas and measure the largest lichen thallus on each one. However, how to express the data in terms of a probability distribution function and link it to an age of last exposure of the rock surface are still a matter of debate. Here, we propose a novel approach to the modeling of lichen growth curves by treating lichen growth as a continuous-time Markov process with a time-varying rate and additive Brownian noise. Given the growth rates, the probability distribution of the lichen population at any time can then be obtained by solving the Fokker–Planck equation. This method is illustrated using a dataset from the Huashan area of eastern China, which consists of measurements of the largest thalli on 12 rock surfaces of known age. We first build up the probability distribution of the lichen population for each rock surface based on extreme value theory and then use these to optimize the growth curve by minimizing the Jensen–Shannon divergence. A new method is also proposed to use the growth curve to map a sample of size data from an undated rock surface to the calendar age domain so as to yield a fully probabilistic estimate of the exposure age of the undated rock surface rather than a point estimate.  相似文献   

ERA-Interim和CMFD气象驱动数据在新疆额尔齐斯河流域的适用性评价     

高黎明  张乐乐  沈永平  张耀南  张伟 《冰川冻土》2022,44(1):179-187

气象驱动数据质量是影响流域水文过程模拟精度的一个重要因素。基于新疆额尔齐斯河流域及周边区域8个气象站记录的数据,对ERA-Interim再分析资料和中国区域地面气象要素驱动数据集（CMFD）在流域的适用性进行了评价,并对比了ERA-Interim和CMFD气象要素年均值在流域的空间分布。结果表明：ERA-Interim和CMFD记录气温、相对湿度、向下短波辐射和向下长波辐射数据与观测数据具有较高的一致性,但降水和风速数据与观测数据的一致性比较差。小时尺度上ERA-Interim记录的气温、相对湿度、降水量、向下短波辐射准确度略高于CMFD数据,而日尺度上CMFD记录的所有气象要素的准确度均高于ERA-Interim数据,结合Noah-MP模型的模拟结果,认为CMFD数据在新疆额尔齐斯河流域的适用性整体优于ERA-Interim数据。从两种驱动数据获取的流域气象要素空间分布来看,ERA-Interim和CMFD获取的年平均气温、风速、相对湿度、降水量、向下长波辐射在流域空间具有高度一致性,但向下短波辐射空间分布差别较大。  相似文献   

三极地区雪冰中黑碳研究进展     

马翔宇  李传金 《冰川冻土》2021,43(1):92-106

黑碳被认为是除温室气体外对气候变暖贡献最大的辐射强迫因子。三极（北极、南极和青藏高原）地区是全球雪冰分布最集中的区域,沉积至雪冰中的黑碳可反映人类活动的历史变化,并可能导致反照率降低而影响物质能量平衡。通过系统回顾三极地区雪冰黑碳的研究方法、空间分布、时间变化及其造成的辐射强迫,得到的研究结果表明： 由于处于不同的地理位置和环境条件,三极雪冰中黑碳的时空分布及其辐射强迫差异较大,其中青藏高原是浓度和辐射强迫影响最大的地区,也是对水资源和生态安全潜在影响最严重的区域。三极地区雪冰中保存着长时间序列的黑碳沉积记录,是研究自然变率、人类活动影响黑碳沉积历史（如北极与青藏高原冰芯中记录的工业革命以来黑碳沉积的快速上升）的良好介质,同样为模型预测未来变化提供了数据支持。三极地区是全球变化的指示器和放大器,在全球变暖不断加剧的背景下,黑碳必然会在未来的三极地区气候演变中扮演更为重要的角色。  相似文献   

GRAPES全球集合预报系统模式扰动随机动能补偿方案初步探究   总被引：10，自引：5，他引：5       下载免费PDF全文

彭飞  李晓莉  陈静  李红祺 《气象学报》2019,77(2):180-195

为了体现次网格尺度能量升尺度转换过程中存在的不确定性, 文中将随机动能补偿(Stochastic Kinetic Energy Backscatter, SKEB)方案应用于GRAPES(Global/Regional Assimilation and Prediction System)全球集合预报系统(GRAPES-GEPS), 以更好地表征模式误差并且增大集合离散度。使用的SKEB方案基于具有一定时、空相关特征的随机型以及由数值扩散导致的局地动能耗散率来构造随机流函数强迫。并根据流函数与水平风速旋转分量的关系, 将SKEB方案中的流函数强迫转化为适用于GRAPES全球模式的水平风速扰动。结果表明, SKEB方案的使用一方面能够提高GRAPES对大气动能谱的模拟能力; 另一方面能够改善GRAPES-GEPS的集合离散度与集合平均误差的关系, 增加了集合离散度, 并在一定程度上减小了集合平均误差, 尤其是在热带地区这种改进更为显著。而且该方案使得热带地区连续分级概率评分(CRPS评分)显著减小。就降水预报而言, 从Brier评分与相对作用特征面积(AROC, Area under the Relative Operating Characteristics)的结果来看, SKEB方案有助于改善中国地区小雨[0.1 mm, 10 mm)、中雨[10 mm, 25 mm)与大雨[25 mm, 50 mm)量级降水的概率预报技巧, 而对暴雨[50 mm, ∞)量级降水预报技巧影响很小(24 h降水量)。总体上, 模式扰动随机动能补偿方案提高了GRAPES-GEPS的概率预报技巧。   相似文献   

Efficient stochastic analysis of unsaturated slopes subjected to various rainfall intensities and patterns     

《地学前缘(英文版)》2023,14(1):101490

Rainfall infiltration poses a disastrous threat to the slope stability in many regions around the world. This paper proposes an extreme gradient boosting (XGBoost)-based stochastic analysis framework to estimate the rainfall-induced slope failure probability. An unsaturated slope under rainfall infiltration in spatially varying soils is selected in this study to investigate the influences of the spatial variability of soil properties (including effective cohesion c′, effective friction angle φ′ and saturated hydraulic conductivity k_s), as well as rainfall intensity and rainfall pattern on the slope failure probability. Results show that the proposed framework in this study is capable of computing the failure probability with accuracy and high efficiency. The spatial variability of k_s cannot be overlooked in the reliability analysis. Otherwise, the rainfall-induced slope failure probability will be underestimated. It is found that the rainfall intensity and rainfall pattern have significant effect on the probability of failure. Moreover, the failure probabilities under various rainfall intensities and patterns can be easily obtained with the aid of the proposed framework, which can provide timely guidance for the landslide emergency management departments.  相似文献   

Geophysical flows under location uncertainty,Part I Random transport and general models     

V. Resseguier  E. Mémin  B. Chapron 《地球物理与天体物理流体动力学》2017,111(3):149-176

A stochastic flow representation is considered with the Eulerian velocity decomposed between a smooth large scale component and a rough small-scale turbulent component. The latter is specified as a random field uncorrelated in time. Subsequently, the material derivative is modified and leads to a stochastic version of the material derivative to include a drift correction, an inhomogeneous and anisotropic diffusion, and a multiplicative noise. As derived, this stochastic transport exhibits a remarkable energy conservation property for any realizations. As demonstrated, this pivotal operator further provides elegant means to derive stochastic formulations of classical representations of geophysical flow dynamics.  相似文献   

Precessional forcing of tropical vegetation carbon storage     

William D. Gosling  Philip B. Holden 《第四纪科学杂志》2011,26(5):463-467

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Geophysical flows under location uncertainty,Part II Quasi-geostrophy and efficient ensemble spreading     

V. Resseguier  E. Mémin  B. Chapron 《地球物理与天体物理流体动力学》2017,111(3):177-208

Models under location uncertainty are derived assuming that a component of the velocity is uncorrelated in time. The material derivative is accordingly modified to include an advection correction, inhomogeneous and anisotropic diffusion terms and a multiplicative noise contribution. In this paper, simplified geophysical dynamics are derived from a Boussinesq model under location uncertainty. Invoking usual scaling approximations and a moderate influence of the subgrid terms, stochastic formulations are obtained for the stratified Quasi-Geostrophy and the Surface Quasi-Geostrophy models. Based on numerical simulations, benefits of the proposed stochastic formalism are demonstrated. A single realization of models under location uncertainty can restore small-scale structures. An ensemble of realizations further helps to assess model error prediction and outperforms perturbed deterministic models by one order of magnitude. Such a high uncertainty quantification skill is of primary interests for assimilation ensemble methods. MATLAB^® code examples are available online.  相似文献   

1D elastic full‐waveform inversion and uncertainty estimation by means of a hybrid genetic algorithm–Gibbs sampler approach          下载免费PDF全文

Mattia Aleardi  Alfredo Mazzotti 《Geophysical Prospecting》2017,65(1):64-85

Stochastic optimization methods, such as genetic algorithms, search for the global minimum of the misfit function within a given parameter range and do not require any calculation of the gradients of the misfit surfaces. More importantly, these methods collect a series of models and associated likelihoods that can be used to estimate the posterior probability distribution. However, because genetic algorithms are not a Markov chain Monte Carlo method, the direct use of the genetic‐algorithm‐sampled models and their associated likelihoods produce a biased estimation of the posterior probability distribution. In contrast, Markov chain Monte Carlo methods, such as the Metropolis–Hastings and Gibbs sampler, provide accurate posterior probability distributions but at considerable computational cost. In this paper, we use a hybrid method that combines the speed of a genetic algorithm to find an optimal solution and the accuracy of a Gibbs sampler to obtain a reliable estimation of the posterior probability distributions. First, we test this method on an analytical function and show that the genetic algorithm method cannot recover the true probability distributions and that it tends to underestimate the true uncertainties. Conversely, combining the genetic algorithm optimization with a Gibbs sampler step enables us to recover the true posterior probability distributions. Then, we demonstrate the applicability of this hybrid method by performing one‐dimensional elastic full‐waveform inversions on synthetic and field data. We also discuss how an appropriate genetic algorithm implementation is essential to attenuate the “genetic drift” effect and to maximize the exploration of the model space. In fact, a wide and efficient exploration of the model space is important not only to avoid entrapment in local minima during the genetic algorithm optimization but also to ensure a reliable estimation of the posterior probability distributions in the subsequent Gibbs sampler step.  相似文献   

On the criteria for the initiation of motion in tidal inlets,deterministic and stochastic approaches     

Morteza Kolahdoozan  Hanifeh Imanian  Roger A. Falconer 《Coastal Engineering》2011

Details are given herein of the refinement and application of a three-dimensional layer integrated numerical model to predict morphological changes in tidal basins. The solution of governing differential equations, which consist of the conservation of mass and momentum for the hydrodynamics, the transport equation for the suspended sediment fluxes and the sediment mass conservation equation for the bed level changes are carried out by the use of Alternating Direction Implicit (ADI) Finite Difference Method (FDM). The model includes different criteria for the initiation of motion namely Shields (1936, Application of Similarity Principles and Turbulence Research to Bed load Movement, Hydrodynamics Laboratory, California Institute of Technology, Pub. No. 167), Kolahdoozan (1999, Numerical Modelling of Geomorphological Processes in Estuarine Waters, PhD Thesis, Department of Civil and Environmental Engineering, University of Bradford, Bradford, UK, 288) and Zanke (2003, On the Influence of Turbulence on the Initiation of Sediment Motion, International Journal of Sediment Research, 18(1), 17–31), to compare different aspects of flow conditions. As the flow is highly turbulent with the random nature of its components, many researchers have tried to express sediment transport processes by using stochastic approaches. In the current study both deterministic and stochastic methods are included in the numerical model to evaluate their accuracy and efficiency. To validate the numerical model results, laboratory measurements are used, with these being obtained from an earlier experimental program undertaken by the authors. Results of a short term bed level changes in a laboratory model harbor are included for the model verification purposes. Comparisons are undertaken using different criteria for the initiation of motion, with the results highlighting that the unsteadiness in the flow parameters included in the numerical model has a major effect on the bed level changes inside the harbor, in compare with the turbulence structure of the flow. The model is then applied to a real case study of the Humber Estuary, located in the UK, with comparisons being undertaken for different criteria for the initiation of motion, using both deterministic and stochastic approaches for the long term bed level predictions.  相似文献